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ABSTRACT 

The  state-space  technique  Is  used  to  formulate  the  transient 
electromagnetic  problem  of  excitation  of  a rectangular  cavity  through 
an  aperture.  Depending  on  whether  the  E-  or  the  H-formulatlon  Is  used, 
the  effect  of  the  aperture  can  be  accounted  for  by  a magnetic  current 
in  Invoking  the  induction  theorem  or  by  an  electric  current  in  in- 
voking the  equivalence  theorem.  The  governing  second-order  differen- 
tial equation  Is  converted  Into  a set  of  first-order  state  equations 
by  defining  three  new  state  variables  In  addition  to  an  appropriate 
vector  potential.  These  state  equations  are  solved  by  the  method  of 
moments.  Two  cases  are  considered:  parallel  polarization  and  perpen- 
dicular polarization.  Some  significant  singularities  for  the 
parallel-polarization  case  are  found  and  the  electric  intensities  as 
functions  of  time  at  two  locations  In  the  cavity  are  computed  for  a 
step  excitation.  An  alternative  approach  of  expanding  the  cavity 
fields  In  terms  of  TM  and  TE  cavity  modes  Is  developed  to  reduce  the 
dimension  of  the  system  matrix.  Numerical  difficulties  are  encountered 
and  possible  contributing  causes  arc  discussed. 
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I . INTRODUCTION 
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Transient  electromagnetic  problems  arise  in  the  consideration  of 
shielding  effectiveness,  electromagnetic  compatibility,  radar  target 
Identification,  remote  environmental  sensing,  and  electromagnetic  effects 
dttt  to  X rays  and  gamma  rays  generated  by  a nuclear  blast.  The  time- 
harmonic  behavior  of  an  antenna  or  a scatterer  can  be  obtained  first  as  a 
function  of  frequency  by  solving  the  governing  Integral  equations  with  the 
given  boundary  conditions.  The  transient  response  can  then  be  determined 
by  performing  an  Inverse  Fourier  transformation.  Numerical  methods  are 
used  in  the  frequency-domain  solution  for  many  different  frequencies.  The 
brute-force  numerical  Fourier  inversion  is  generally  inefficient  and  con- 
vergence problems  arise  in  superposing  the  steady-state  solutions.  A few 
simple  situations,  such  as  transient  scattering  from  wire  antennas  and  con- 
ducting cylinders  Il]-[3],  have  been  analyzed  directly  in  the  time  domain. 

i 

Basls-fimctlon  expansions  and  inner  products  over  both  space  and  time  are  I 

i 

required.  Space-time  Integro-dlfferential  equations  are  encountered  and 
the  numerical  representations  of  the  derivatives  lead  to  very  complicated 
procedures . 

A variation  of  the  time-domain  approach  makes  use  of  a a Hall^n-type  ^ 

integral-equation  formulation  [4].  As  a consequence  of  the  absence  of  space- 
time  derivatives  imder  the  Integrals,  the  numerical  process  is  less  compli- 
cated. However,  because  of  retarded  time,  a tedious  step-by-step  iterative 
procedure  has  to  be  used  in  conjunction  with  the  proper  families  of  trajec- 
tories known  ao  characteristic  curves  in  order  to  determine  the  homogeneous 
solutions.  All  of  the  methods  mentioned  above  are  tedious  and  do  not  afford 

3 


• physical  Insight  in  the  solution  of  a transient  problem.  Moreover,  a 
complete  recalculation  would  be  necessary  under  any  change  In  the  wave 
shape,  polarization,  or  angle  of  Incidence  of  the  source  of  excitation. 


More  recently  a singularity-expansion  method  has  been  used  to  deter- 
mine the  transient  response  of  scatterers  of  simple  geometries  [5]-[ll]. 
Responses  to  tr£inslent  excitations  are  expressed  In  terms  of  exterior 
natural  frequencies,  modes,  and  coupling  coefficients,  and  Induced  currents 
are  represented  by  a series  of  damped  sinusoidal  functions.  This  method  has 
the  advantages  of  providing  a physical  sight  to  the  radiation  or  scattering 
problem  and  of  allowing  the  response  to  be  determined  for  a change  In  source 
parameters  without  a complete  recalculation.  In  the  evaluation  of  the  natural- 
mode and  coupling  vectors.  It  Is  necessary  to  know  the  nature  of  the  singu- 
larities Involved  in  the  inverse  Laplace  transform.  Knowledge  in  this  respect 
is  not  yet  secure.  The  transcendental  nature  of  its  system  Impedance  matrix 
results  In  an  Infinite  number  of  complex  poles  whose  locations  must  be  numeri- 
cally searched.  It  Is  generally  agreed  that  the  singularity-expansion  method 
Is  not  very  satisfactory  for  evaluating  early-tlme  responses. 

The  method  of  characteristic  modes  has  been  used  for  determining  the 
steady-state  response  of  conducting  bodies  [12].  A striking  similarity  appears 
to  exist  between,  this  method  and  the  singularity-expansion  formulation  if  the 
method  of  moments  [13]  using  a common  spatial  basis  Is  applied  to  both  cases; 
but  no  formal  relations  have  yet  been  established.  The  natural  modes  from 
the  singularity  expansion  are  not  orthogonal.  The  characteristic  modes  are 
orthogonal,  but  they  vary  with  the  source  frequency  and  it  would  be  necessary 
to  compute  the  characteristic  modes  for  all  frequencies  before  a Fourier  Inver- 
sion could  be  effected  to  determine  the  transient  behavior. 
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The  problem  of  transient  field  behavior  Inside  a conducting  cavity 
due  to  excitation  through  an  aperture  by  an  Incident  electromagnetic  pulse 
(EMP)  is  particularly  difficult  because  of  the  reflections  from  the  cavity 
vails  and  the  coupling  between  the  Interior  and  exterior  fields  at  the  aper- 
ture. Past  Investigations  on  EMP  excitation  of  cavity-backed  apertures  have 
largely  dealt  with  small  openings  and  have  neglected  the  effect  of  cavity 
reflections  on  the  aperture  field  distribution.  For  small  openings  the 
quasl-statlc  method  Is  used  to  determine  the  fictitious  magnetic  current 
and  charge  distributions  In  the  aperture.  Equivalent  electric  and  magnetic 
dipoles  are  defined,  and  their  radiated  fields  determined  with  the  aid  of 
scalar  and  vector  potentials  [1A]-[19].  The  fields  In  the  cavity  are  custom- 
arily expanded  In  terms  of  unperturbed  normal  modes.  The  quasl-statlc  approxi- 
mation cannot  be  applied  when  the  aperture  is  not  small  and  when  early-time 
responses  are  Important.  In  neglecting  the  effect  of  the  reflections  from 
caAd-ty  walls  on  the  aperture  field  distribution,  one  essentially  treats  the 
external  and  Internal  portions  of  the  problem  separately.  Since  cavity  dimen- 
sions obviously  play  an  Important  part  in  the  total  problem,  this  approach  may 
result  In  significant  errors. 

In  this  report  we  will  avoid  the  quasl-statlc  approximation  and  solve 
the  internal  and  external  portions  of  the  problem  simultaneously.  New 
variables  (state  variables)  will  be  Introduced  to  convert  the  governing 
second-order  differential  equation  into  a set  of  first-order  equations  which 
correspond  to  normalized  state  equations.  The  field  within  the  cavity  will 
be  expanded  In  terms  of  suitably  chosen  subsectional  expansion  functions  with 
variable  coefficients  and  the  field  outside  the  cavity  expressed  as  a super- 
position of  plane-wave  fields.  The  cavity  and  the  external  fields  are  matched 
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at  Che  aperture  where  a fictitious  equivalent  current  exists.  A combined 
field  expression  containing  the  unknown  expansion  coefficients  Is  obtained. 

To  determine  these  coefficients  the  method  of  moments  [13]  Is  used  to  con- 
vert the  first-order  equations  Into  matrix  equations.  It  will  be  shown  that 
the  typical  coefficient  matrix  can  be  expressed  In  a form  for  which  the 
singularity-expansion  method  [5]  can  be  used  to  advantage. 

The  general  procedure  of  solution  for  cavity-backed  aperture  problems 
Is  outlined  first.  The  theoretical  formulation  for  the  transient  excitation 
of  a rectangular  cavity  with  a slot  aperture  Is  then  given  for  both  parallel 
and  perpendicular  polarizations.  Some  numerical  results  are  Included. 
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II . SOLUTIOW  PROCEDURE 
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We  consider  the  problem  of  a slot  aperture  In  an  Infinite  conducting 
plane  backed  by  a rectangular  conducting  box,  as  shown  in  Fig.  1.  An  inci- 
dent transient  electromagnetic  wave  (E^,  H^)  impinges  normally  on  the  plane 
an4  the  aperture.  The  problem  is  to  determine  the  scattered  field  in  the 
y > 0 region  and  the  field  penetrated  through  the  aperture  into  the  conduct- 
ing cavity. 

The  induction  theorem  [20],  [21]  can  be  invoked  for  the  solution  of 
this  problem.  Figure  2(a)  represents  a simplified  2-dlmenslonal  view  of 
the  original  problem.  (E  , H ) and  (E  , H.)  are,  respectively,  the  cavity 
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field  and  the  external  scattered  field.  In  order  to  determine  these  unknown 
fields,  we  consider  the  case  when  the  aperture  is  covered  by  a conductor.  The 
entire  region  to  the  left  of  the  infinite  plane  will  have  a null  field  and, 
according  to  the  induction  theorem,  a magnetic  current  on  the  right  surface 
of  the  conducting  plane  will  support  a different  scattered  field  (E®,  H°),  as 
shown  in  Fig.  2(b),  where 

M = E°  X ft 
o s 

- ft  X 

- ^ X . (1) 

For  a normally  incident  plane  wave  (E^,  H^) , the  scattered  field  (E°,  5°) 

from  an  infinite  conducting  plane  without  an  aperture  is  easily  determined. 

The  null  field  to  the  left  of  the  plane  will  be  maintained  if  the  plane  is 

removed  and  a magnetic  current  2M  exists  in  its  place  which  will  result  in 

o 

a field  (E^  + E°,  + H*)  in  the  y > 0 region,  as  shown  in  Fig.  2(c). 
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Subtracting  the  fields  In  Fig.  2(c)  from  those  in  Fig.  2(a),  we 
obtain  tb«.  problem  in  Fig,  2(d).  The  magnetic  current  M in  the  aperture 

Is 

M - - 2M^  - - 29  X , (2) 

which  supports  the  field  (E  , H ) inside  the  cavity  and  a field  (E  - 1°, 

C C 8 S 

Hg”  H^)  to  the  right  of  the  infinite  plane.  We  note  that  the  region  in 

which  the  difference  field  (E  - E°,  H - H°)  exists  is  source-free  and 

s s s 8 

that  the  tangential  component  of  the  electric  field  is  required  to  vanish 
on  conducting  walls. 

For  the  problem  in  Fig.  2(d),  we  start  from  the  two  Maxwell's  curl 
equations 


7 X E 


- % 37  “ “ 


V X H = e 


M. 

O dt 


(3) 

(4) 


Taking  the  Laplace  transform  of  Eqs.  (3)  and  (4),  we  obtain 


7xe*-psH-M  (5) 

o 

7 X S = 0^8  I , (6) 

where  a tilde  (-)  over  a quantity  denotes  the  Laplace  transform  of  that 
quantity. 

Let  P be  the  Laplace  transform  of  an  electric  vector  potential  F 
such  that 

E = - 7 X F . (7) 

Combining  Eqs.  (5)  to  (7)  and  using  the  Lorentz  gauge,  we  have  an  in- 
homogeneous Helmholtz  equation: 
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(8) 


-2~  2 — - 
V F - w e 8 F - - M . 
o o 

31^  ^ 

Solution  of  Eq.  (8)  for  F will  give  R from  Eq.  (7)  and  H from 

H - — V X 7 X f 


I 

(9) 
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f 

r 

I 

( 
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In  regions  where  M Is  zero. 

Specialization  of  these  general  formulas  will  depend  on  the  polari- 
zation of  the  Incident  wave;  but  as  soon  as  E^  Is  known,  Eqs.  (7)  and  (9) 
can  be  expanded  Into  component  equations  and  the  source  term  R In  Eq.  (8) 
can  be  found  from  Eq.  (2).  A set  of  new  variables  can  then  be  defined 
which  will  convert  the  second-order  differential  equation  (8)  Into  a set  of 
first-order  equations,  and  these  equations  are  Laplace-transformed  normalized 
state  equations  In  the  new  state  variables.  Solution  of  the  transformed 
state  equations  Involves  four  steps.  First,  the  space  inside  the  cavity  Is 
divided  Into  subsections  and  suitable  expansion  functions  are  chosen  over 
the  subsections.  The  elements  of  the  unknown  state-variable  vector  are  then 
expressed  In  terms  of  the  expansion  functions  within  the  cavity.  Second,  the 
field  in  the  y > 0 region  Is  expressed  as  a superposition  of  plane  waves. 
Third,  the  cavity  and  the  half-space  fields  are  matched  at  the  aperture. 
Fourth,  inner  products  are  taken  so  that  the  matrix  equations  for  the  un- 
known expansion  coefficients  are  obtained.  These  steps  are  outlined  sepa- 
rately for  the  cases  of  parallel  and  perpendicular  excitations  in  following 
sections. 


III.  RECTANGULAR  CAVITY  WITH  SLOT  APERTURE 


PARALLEL  POLARIZATION 


In  this  section  we  consider  the  case  of  an  Incident  plane  wave  with 
the  electric  field  polarized  In  a direction  parallel  to  the  slot.  Referring 
to  Fig.  1,  we  have 

- 2 (10) 

z 

and  the  Laplace  transform  of  Eq.  (2)  becomes 

M - - 2R  E^  - R M (11) 

z X 

which  has  only  an  x-component.  The  x-component  of  Eq.  (8)  Is  then 


V^F  - y e 8^  F “ - M 6(y) 

X O O X X f 


(12) 


where  6(y)  is  a Dirac  delta  fimction.  From  Eqs.  (7)  and  (9),  we  have 


E = 0 

X 


~ F 
3z  X 


E = F 
z 3y  X 

ii  = - — [-^  + -^]F 

X y s \ 2 . 2 X 

o 3v  3z 


1 3 


y y^s  3x3y  x 


H 


1 3 


z y s 3x3z  X 
o 


F . 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


The  second-order  differential  equation  (12)  can  be  represented  as  a 
set  of  first-order  equations  by  defining  new  quantities  u,  v,  and  w such 
that 


10 


(19) 


3 ** 

3^  F^(r,8)  - 8 u(r,8) 


and 


* y'-»> 


8 v(r,8) 


8 w(r,8)  , 


(20) 

(21) 


where  r Is  the  space  variable.  We  have,  from  Eq.  (12), 


— Q(r,s)  + ^ ^(r,s)  + ^0(r,8>  - - Fj^(r»8)  + ^ \ • (22) 

Comparing  Eqs.  (21)  and  (20)  with  Eqs.  (14)  and  (13)  respectively,  we  see 
that 


and 


E - 8 « 

y 


E “ - 8 V 
z 


(23) 

(24) 


The  Introduction  of  u,  v,  and  w and  the  use  of  the  first-order  equations 
%d.ll  result  In  significantly  faster  convergence  In  the  numerical  solution. 

The  first-order  equations  (19)  to  (22)  can  be  written  In  a succinct 
form  by  defining  the  following  operators  and  column  matrices: 


_3 

9x 


_9 

3y 


_3_ 

9z 


2 

3x 


0 -i. 


3y 


3z 


- 0 


- 0 


0 

0 


0 

0 


(25) 
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p - 


P e 
o o 


0 

1 


f(r,8)  - 


and 


Equations  (19)  to  (22)  become 


u- 

V 

w 


- M 

S X 


0 

0 

0 


0 

0 


0 

0 


L f(r,s)  = s P f(r,s)  + e (r,s)  . 
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(26) 


(27) 


(28) 


(29) 


Note  that  the  inverse  Laplace  transform  of  Eq.  (29)  is  a set  of  normalized 
state  equations  in  the  four  state  variables  F^,  u,  v,  and  w. 

In  order  to  solve  Eq.  (29)  by  the  method  of  moments  we  subdivide 
the  space  within  the  cavity  in  the  x,v>  and  z directions  and  choose  expan- 
•Ion  function.  (t) , «.d  ^^(c)  ouer 

the  siibsectlons.  The  expansion  functions  must  satisfy  the  required  boundary 
conditions.  For  convenience,  we  define  the  following  column  vectors: 
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*'x(i,j,k)^’'^ 

0 

0 

0 

• 'u.J.k)<">- 

“(I.J.k)'"’ 

0 

0 

0 

0 

0 

0 

'(I.J.k)'"’  * 

0 

'(I.J.k)'"’  • 

\l.J.k)^"> 

0 

0 

/(I.J.k) 

In  view  of  Eq.  (27),  we  can  then  write  the  expanded  form  of  f(r,s)  Inside 
the  cavity  as 


(i.J.k) 


(s)  f 


(l,J,k) 


(r)}. 


Note  that  the  expansion  functions  F^,  u,  v,  and  w are  functions  of  position 
only  and  that  the  Inverse  transformation  of  5,  and  J will  yield  the 

time-varying  expansion  coefficients. 

The  field  in  the  half-space  y 5^  0 Is  expressed  as  a superposition  of 
plane  waves  and  the  Internal  and  external  fields  are  matched  at  the  aperture 
y ■ 0.  A combined  field  expression  for  l(r,s)  can  be  obtained  which  holds 
Inside  the  cavity.  In  the  y > 0 half-space,  as  well  as  In  the  slot.  This 
Is  then  substituted  In  Eq.  (31)  which,  after  Inner  products  with  the  expan- 
sion functions  In  Eq.  (30)  have  been  taken,  leads  to  the  following  matrix 
eqtiatlon: 


(31) 
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where  fi  i **’®  expansion  coefficients  for  the  magnetic  current  M , 

Xvi  ,ny,K  ) X 

C(s)  Is  a column  matrix  and  I 's,  p 's  and  q's  are  themselves  matrices 

mn  mn 

arising  from  the  inner  products,  m and  n are  Indices  locating  the  position 
of  a particular  siibsection  over  which  an  Inner  product  Is  taken.  The  expres- 
sions for  C(s),  Pmn'®  very  complicated.  They  have  been 

given  In  a previous  report  [22],  and  will  not  be  repeated  here. 

The  unknown  coefficient  matrices  {a(s)},  {B(s)},  {y(s)}  and  (8(8)}  can 
be  solved  from  Eq.  (32).  An  outline  of  the  method  of  solution  will  be  given 
In  Section  V.  Determination  of  these  coefficient  matrices  enables  the  calcu- 
lation of  f(r,s)  In  Eq.  (27)  from  Eq.  (31),  which.  In  turn,  leads  to  the  field 
inside  the  cavity. 


IV.  RECTANGULAR  CAVITY  WITH  SLOT  APERTURE 


- PERPENDICULAR  POLARIZATION 


'When  the  electric  field  of  the  incident  plane  wave  is  polarized  in 
• direction  normal  to  the  slot  shown  in  Fig.  1;  that  is,  if 

, (33) 

z 

It  is  store  contrenient  to  use  an  H-formulatlon.  Instead  of  invoking  the 
Induction  theorem,  we  apply  the  equivalence  theorem  [20],  [21],  The  H- 
fleld  discontinuity  at  the  aperture  in  Fig.  2(d)  is  supported  by  an  elec- 


tric current 


J - - 2<l  X 


- 2R  H “ R J . 

z X 


To  solve  this  problem,  we  start  with  the  following  Maxwell's  equations 


which  transform  to 


V X E 


Vx|--vsH 

o 


Vxh-csE  + J. 
o 


Let  A be  the  Laplace  transform  of  a magnetic  ■vector  potential  A such  that 


H - V X A 


We  ha've,  from  Eqs.  (37)  to*(39). 
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V^A 

X 

- p e s^A  ■ -”J  fi(y) 

(40) 

H - 

X 

0 

(41) 

H - 

y 

^ A 

3z  X 

(42) 

H - 

z 

-f- A 

3y  X 

. »2  «2 

(43) 

E - 

X 

L.  r^_  + -i— ]A 

V W * 

(44) 

E - 

y 

e 8 dxdy  X 

0 

(45) 

E - 
z 

_L-  ^ 

e 8 3x3z  X 

(46) 

The  second-order  differential  equation  (40)  can  be  converted  Into 
a set  of  first-order  equations  by  defining  new  quantities  Q,  y,  and  O 
similar  to  those  In  Eqs.  (19)  to  (21)  with  replaced  by  A^.  Instead  of 
Eqs.  (22)  to  (24),  we  now  have 


H - - 8 « (48) 

y 

it  - 8 ^ . (49) 

Z 

An  operator  equation  (29)  which  represents  Laplace  transformed  normalized 
state  equations  In  the  four  state  variables  A^,  u,  v,  and  w Is  obtained 
when  operators  L and  P and  column  matrices  f(r,8)  and  ^gCs)  are  defined  as 
in  Eqs,  (25)  to  (28)  respectively  with  replaced  by  A^  and  replaced 
by  J^.  This  operator  equation  can  be  solved  by  the  method  of  moments  as 
before  by  choosing  suitable  expansion  functions  j j (f) » 

V,.  . !,%(*')»  A over  cavity  subsections.  We  can  write  the 
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'expanded  form  of  l(r,8)  Inside  the  cavity  as  [22] 


’d.j.k)'-'  ‘d.j.k) 


i.j.k  ' "(i.j.k) 


where  the  column  matrices  ^(i^j » ^(i  j.k)^*^^*  ^(1  j *” 

the  same  as  those  defined  In  Eq.  (30), 


0 

0 

0 


(51) 


^d  S,  Y,  and  5 are  the  Laplace  transforms  of  the  expansion  coefficients. 
Analogously  to  the  parallel-polarization  case  treated  in  Section  III, 


we  can  expand  H in  the  half-space  y ^ 0 as  a superposition  of  plane  waves: 
Z 


-j (k^xfk  y+k  z) 

y ® dk  dk  , 

X z 


(52) 


where 


Jk^  - /(s/c)^  + k;^  + k^ 
J X z 


(53) 


H, 


and  the  new  quantity  g (,  ,k  ) can  be  determined  from  the  boundary  con- 

X Z 


dltlon  at  the  slot  [22].  There  Is  a discontinuity  In  H at  the  slot  on 

z 


account  of  the  equivalent  electric  current  which  can  be  expanded  as 


ixW.o  .k)(»)  ''(l.n  “•  •> 

y y 


(54) 
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A system  of  three  equations  are  obtained  by  matching  Q,  and  it  at 
the  slot.  A combined  field  expression  which  holds  Inside  the  cavity. 

In  the  y > 0 half-space,  as  well  as  In  the  slot  can  be  written.  This 
expression  Is  extremely  complicated.  Suffice  It  to  say  that  when  It 
is  substituted  In  Eq.  (50)  and  Inner  products  with  the  expansion  func- 
tions are  taken,  we  obtain  a matrix  equation  similar  to  Eq.  (32),  from 
ifhlch  the  unknown  coefficient  matrices  {a(s)},  {B(s)},  {y(8))»  and  {^(s)} 
can  be  solved. 


V.  SOME  NUMERICAL  RESULTS 


In  this  section  we  outline  the  procedure  for  determining  the  unknown 

coefficient  matrices  {a(s)},  {0(s)},  (yCs)},  and  {6(s)}  from  Eq.  C32) . Al- 

^ j /nUV  -vu  .w  -vw  Fv  uv  vF  vu  w 

though  the  component  matrices  (£  , 1 , I . 1 . p , p , p , p , p , 

mn  mn  mn  mn  mn  mn  mn  mn  mn 

w 

and  p ) representing  the  coupling  between  cavity  and  external  fields  appear 
hl{^ly  complex,  they  are  relatively  sparse.  Typically  an  equation  of  the 
following  form  is  obtained  from  Eq.  (32): 


[Z(8)l  {a(s)}  = [H(s)]  ^ 

where  [Z(s)]  and  [H(s)]  are  square  matrices  containing  I 's,  p 's,  and 

mn  mn 

q *8  In  Eq.  (32),  and  {K(s)}  Is  In  general  not  the  same  as  {C(s)}.  From 

Eq.  (55)  we  have 


{5(8)}  = [Z(s)]“^  ([H(s)]  „ .x(s)}  + {K(8)})  . (56) 

s x(i,ny,k; 

Similar  expressions  are  obtained  for  {6(s)},  {•y(s)},  and  {{(s)}  . 

Let  s^  be  the  zeros  of  |z(s)|  or  the  roots  of  the  equation 

det[Z(s))  = 0 . (57) 

In  circuit-theory  terminology,  [Z(s)]  corresponds  to  the  system  Impedance 

matrix  and  s are  the  natural  frequencies.  [Z(s)]  ^ can  be  expanded  In  a 
a 

partial-fraction  form  as  follows: 

, [R  ] 

[Z(8)r^  = I » (58) 

a a 

%rhere  the  constant  square  matrix  [R^]  Is  the  system  residue  matrix  at  the 

pole  s . (R  1 can  be  written  as  the  product  of  a natural  mode  vector  {R™} 
' a a a 

and  the  transpose  of  a coupling  vector  {R^}  [5],  [9]: 
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I t 


'“.I  ■ "!>  • 


where  {R™}  is  a solution  of  the  equation 


[Z(8  )]  {R“}  - 0 
a a 


and  {R^}  Is  a solution  of 


(Z(s  {R^}  - 0 . 

a a 


A close  examination  of  the  composition  of  the  matrices  [ZCs)]  and 
[H(8)]  reveals  that  their  poles  coincide  and  therefore  cancel.  We  have, 
from  Eqs.  (56),  (58)  and  (59), 

{r®}{rC}T 

{a(8)}  - I — ((H(s)]  {i  A ^ _(8)}  + {K(8)})  . (62) 

a » y*  ' 

Now  define 

IH(8)]  „ , .(s)}  + {K(s)}  = N(s)  {V^(8)}  , (63) 

8 xCi,ny,Ky  o 

where  {V  (s)}  Is  the  excitation  vector  when  the  incident  wave  is  a pulse, 
o 

We  can  then  write  Eq.  (62)  as 

{5(8)}  - I — N(8)  {v^(8)} 

a “a 


where 


a a 


h (s)  - {r;}’'{V^(b)} 
o a o 


is  called  the  coupling  coefficient  [5].  We  note  that  N(s)  itself  may  con- 
tain poles  in  the  finite  plane,  but  this  fact  does  not  result  in  any  serious 
difficulty. 
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We  are  now  In  a position  to  write  the  expressions  for  the  field 


distributions  within  the  cavity.  Vrom  F,q.  (31), 

Vx,y,z,.)  - 

which,  in  view  of  Eq.  (64),  becomes 


(66) 


Fjj(x.y.z,s)  - I n^(s)  {F^^^j^^^(x,y,z)}  (s  - s^)“^  N(s) 


-1 


I hj|(8)  v^Cx,y,2)  (s  - 8^)"^  N(s) 
a 


In  Eq.  (67), 


(67) 


is  a natural  mode  for  F . In  a similar  manner,  we  will  eet 

X ’ 


fi  (x,y,z,8)  - I f^^(8)  vJ(Xiy,z)  (s  - s^)“^  N(8) 


and 


E (x,y,z,8)  » I n (s)  V ^(x,y,z)  (s  - s ) ^ N(8) 
Z C(  a a 


(68) 


(69) 


(70) 


EE 

where  v ^ and  v * are  the  natural  modes  for  E and  £ respectively, 
a a y z 

For  numerical  computation,  a rectangular  cavity  of  dimensions 
4x6x2  was  chosen  and  the  slot  width  was  a/10  or  0.4.  Considerable  dif- 
ficulties were  experienced  in  the  determination  of  the  singularities 
(natural  frequencies)  s^.  After  having  compared  several  root-finding 
procedures  and  checked  meticulously  our  computer  programs,  we  obtained 
the  locations  of  the  singularities  in  the  upper  part  of  the  left  half- 
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plane  for  the  parallel-polarization  case  with  E-formulatlon  as  shown  In 
Fig.  3.  Two  layers  of  singularities  were  found:  the  first  layer  was  very 
close  to  the  Ja>-axis  and  the  second  layer  had  a positive  slope.  The  singu- 
larities in  the  first  layer  fall  on  a smooth  locus,  and  they  appear  to  be 
Irregularly  spaced.  No  method  exists  which  could  ascertain  the  effect  of 
numerical  noise,  subdivision  scheme,  etc.  It  is  known  that  numerical  re- 
sults are  sensitive  to  round-off -errors  and  may  diverge  with  an  increasing 
sampling  density  [23].  Problems  of  numerical  anomaly  appear  in  even  the 
simplest  cases  [23] , and  precluded  the  determination  of  the  natural  fre- 
quencies of  some  well-defined  geometric  structures  such  as  oblate  sphe- 
roids [2A]. 

In  our  ccalculatlons  eight  subsections  were  used  in  the  y-dlrection 
of  the  cavity  and  six  subsections  in  the  x-dlrection.  The  expansion  (and 
testing)  functions  in  the  y-  and  x-  directions  were,  respectively,  rec- 
tangular and  triangular  functions.  The  natural  modes  of  the  first  several 
singularities  of  the  first  layer  were  determined.  Typical  amplitude  and 
phase  plots  at  two  locations  in  the  cavity  for  the  first  pole, 

Sj^j^  «=  (-0.110  + jl.258)c,  and  for  the  third  pole,  Sj^^  ” (-0.141  + j2.296)c, 
are  given  in  Figs.  4 and  5 respectively.  Assuming  a step  excitation,  the 
normalized  cavity  field  |E  |/E^  Is  shown  in  Fig.  6 as  a function  of  time 
at  the  point  (2,  -5,  1),  which  is  (5/6)th  of  the  way  in  from  the  slot 
aperture.  Figure  7 shows  a similar  plot  at  the  point  (2,  -3,  1).  The 
figures  indicate  that  the  maximum  value  of  jE^I  inside  the  cavity  is  in 
the  order  of  one-thousandth  of  the  incident  field  intensity,  E^,  for 
parallel  polarization  (good  shielding  effectiveness) . Maxima  and  nulls 
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exist  as  t varies.  However,  a physical  interpretation  of  their  spaclngs 
is  difficult  because  of  the  complicated  nature  that  the  penetrated  field 
is  reflected  from  the  cavity  walls. 
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VI.  AN  ALTERNATIVE  APPROACH 


r' 


The  formulations  In  the  previous  sections  required  the  expansion 
of  the  cavity  fields  In  terms  of  the  selected  expansion  functions.  As 
a consequence,  the  three-dimensional  moment  method  leads  to  system 
Impedance  matrices,  [Z(s)],  of  large  dimensions  which  often  tend  to 
be  Ill-conditioned.  With  a view  to  avoiding  numerical  difficulties 
Inherent  with  large  matrices,  an  alternative  approach  Is  Investigated 
In  this  section.  Briefly,  the  fields  inside  the  cavity  are  expressed 
in  terms  of  TM  and  TE  cavity  modes  with  unknown  coefficients.  This 
approach  effectively  changes  the  three-dimensional  moment  method  into 
a two-dimensional  one,  thus  reducing  the  dimensions  of  the  resulting 
matrix.  We  develop  this  alternative  approach  for  the  perpendicular- 
polarization  case  with  the  incident  normal  to  the  slot.  As  shown 
in  Section  IV,  the  H-field  discontinuity  at  the  aperture  is  supported 
by  an  equivalent  current  J given  in  Eq.  (34). 


VI -1.  Cavity  Modes 


Referring  to  Fig.  1,  we  separate  the  fields  Inside  the  cavity 

into  TM  to  z (H  =0)  and  TE  to  z (E  = 0)  modes  and  express  the  field 
z z 

TM  TE 

components  in  terms  of  wave  functions  <{>^  and  (j)  respectively. 

I.  TM  to  z mode 


E 

X 


- I™ 

y 3x3z 


i 3^  iTM 
5^  3y3z  ^ 


(71) 

(72) 


jj  . 1 . 32  rTM 

^ c 


(73) 


3 rTM 


A 

X 3y 


5 L ;TM 

“y  3x  * 


H = 0 
z 


(75) 

(76) 


where 


5^  « se. 


•"TM 

The  transformed  wave  function  (p  can  be  written  as 


(77) 


tTM 


where 


y A (s)  sin  k X cos  k z sin  k (y+b) 
mn  X z y ■' 

m,n  ^ 


k =21 
X a 

k “ — 
z c 


(78) 

(79) 

(80) 


and 


c 

o 

rTM 


(81) 


Note  from  Eq.  (78)  that  $ = 0 at  y = - b.  in  8<l*  (81)  Is  used  for 

the  speed  of  light  In  order  to  distinguish  It  from  the  cavity  dimen- 
sion c.  Substituting  Eq.  (78)  In  Eqs.  (71)-(73)  and  (76),  we  have 


fi  ■ 77  y (-1)  k k A (s)  cos  k X sin  k z sin  k (y+b)  (82) 

X 9 X z mn  x z y 

■'  m,n 


E = y (-1)  k k A (s)  sin  k x sin  k z cos  k (y+b)  (83) 

y y m n yzmn  x z y'-'  ' 


^ (-1)  (k^  + ■^)  \n(s)  sin  k^x  cos 
m.n  c^ 

H - 0 
z 


k z sin  k (y+b)  (84) 

z y 


(85) 
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The  wave  function  satisfying  the  boundary  condition  at  y = -b  Is 


Finally  we  combine  the  E and  fi  expressions  for  the  TM  and  TE 

z 

modes  In  Eqs.  (82)-(85)  and  (94) -(97)  to  obtain 


fi,  - I (-l)Ij  l^yCy+b)  (98) 


m,n 


®y  “ I ‘^°®  l^y(y+b) 

m^n  ^ 


(99) 


Sj  “ I (-1)  ~ (fcg  + ^)  ^jjjj(s)  sin  k^x  cos  k^z  sin  k (y+b) 
fflf  n •“  ^ 


and 


o 

2 


• 1 * 

H ■ J (-1)  — (k  + ■^)  8 „(s)  cos  k X sin  k z cos  k (y+b) 

z “ z z z fnn  V «T  ^ 


2 ' z 2'  mn' 
m.n  c^ 


(100) 


(101) 


VI-2«  External  Fields 

The  fields  outside  the  cavity  in  the  y ^ 0 region  can  be  expressed 

as  a superposition  of  plane  waves  with  coefficients  determined  by  the 

tangential  components  of  the  E-field  at  the  aperture.  The  two  tangential 

components,  E and  E , are  to  be  treated  separately. 

X z 

(a)  - 0,  E^  ,«  0 

Let 


y>0 


dk'  E (k’,k')e 

Z X X z 


-j(k'x  + k'y  + k'z) 
X y'  z 


(102) 


where 


jk'  - 

o 


(103) 


and 
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Choose  a wave  function 


* r-;  f 


• -J(k>  + k'y  + k'z) 

dk;  f (k^,  kpe  * y 


(105) 


The  Fourier  transform  of  if  Is 


-Jk'y 

4;  - f(k'  kMe  y 

A Z 


(106) 


The  field  components  In  terms  of  ip  are  [21] 


= jk^ip 


Ey  = - jk^ip 


E = 0 
z 


H 


k'k’ 

X z 7 
4) 


k'k' 

fl 2^iP 

y 


(s/c  )^  + k'^ 


H 


(107) 

(108) 

(109) 

(110) 

(111) 

(112) 


Combining  Eqs.  (102),  (107)  and  (112),  we  obtain 


ft  (x,y,z) 


, f»  j(s^/c^+k'^)  -j(k'x  + k'y  + k'z) 

dk'l  dk:'  - 8,(k:.k:)e ' * y"  ^ 


z sp  k' 

o y 


X X z 


(113) 


It  Is  noted  that  £ (k',  k')  can  be  determined  from  Eq.  (104)  for 
XX  z 

a known  £^(x,o,z)  at  the  aperture. 


(b)  - 0,  0. 


In  this  case  we  write 


y>^0 


dk'  E (k',  k')e 

Z Z X z 


-J(k'x  + k'y  + k'z) 

JN  X y' 


(llA) 


where 


Ez(k^.  k^)  * I dx  I dz  g^(x,o,z) 


j(k^x  + k^z) 


(115) 


By  a similar  procedure,  we  obtain  H (x,y,z)  in  the  y^O  region  as 

Z 


H (x,y,z)  = 
z 


-jk’  k' 

X Z.5 


4x 


•CO  *00 

K dk’  dk’  ( H;t)E  (k’,  k’)e 

J X J z ' sp_  k’/  z'  x’  z' 


-j(k^x  + k^y  + k^z) 


(116) 


o y 


Combining  Eqs.  (113)  and  (116),  we  have 


H (x,y,z) 
z 


^ r-i  f 


k’k’E  (k’,k’)  - (s^/c^+k’^)E  (k’k’) 


dk’  {• 
z 


X Z Z X'  z 


O Z X X z 


m^CO  <*00 


sp  k’ 

o y 

-j(k’x  + k’y  + k’z) 
X y'  z 


(117) 


Vl-3.  Boundary  Conditions  at  Aperture 

The  boundary  conditions  at  the  aperture  are  such  that  and 

g are  continuous  and  that  fi  Is  discontinuous  by  the  equivalent 
X 2 

electric  current  J^.  We  treat  the  three  components  separately 
below. 
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(a)  Boundary  condition  for  E^. 

E Is  continuous  at  the  aperture,  y*0.  Applying  the  method  of 
2 

moments,  we  expand  E at  y*0  over  the  aperture  region  In  terms  of  appro- 

* E^ 

prlate  expansion  functions  f with  unknown  coefficients  B (s) . 

pq  pq 

E 


E 


ly=0 


1 aperture 

ptq  ” 


elsewhere 


(118) 


Equating  Eq.  (118)  to  the  of  the  cavity  field  In  Eq.  (100)  at  y=0. 


we  have 


y (-1)  (k^  + •^)A  (s)  sin  k x cos  k z sin  k b 

^ y z 2 mn  X z y 

m,n  c^ 


I * 

p.,  M P<1 


Now  we  take  the  Inner  product  of  Eq.  (119)  with  testing  functions  and 
the  result  can  be  put  In  the  following  form. 


au)}  = [fi]{g  (s)} 
mn  pq 


(119) 


(120) 


It  Is  important  to  note  that  the  dimension  of  the  matrix  [M]  Is  not  very 

large  because  neither  the  number  of  cavity  modes  (m,n)  nor  the  number  of 

subsections  over  the  aperture  (p,q)  needs  to  be  large.  E (k',  k')  can  be 

2 X 2 

obtained  from  Eq.  (115). 
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(b)  Boundary  condition  for  E^, 


E is  also  continuous  at  the  aperture.  E at  y*0  can  be  expanded 

E * 

in  terms  of  expansion  functions  f 

****  E 

^ aperture 

p.q 


|y-0 


(121) 


elsewhere 


Equating  Eq.  (121)  to  the  E^  of  the  cavity  field  in  Eq.  (98)  at  y-0. 


we  have 


I (-l)[i  ^x*  ®^“  V 

m,n  ^ ^ ^ 


I Snn<®)  C 
mtn  P‘1  P‘1 


(122) 


Inner  product  of  Eq.  (122)  with  testing  functions  yields 


{B  (s)}  = (s)}  + [N,I{g  (s)} 

mn  1 pq  2 pq 


(123) 


(c)  Boundary  condition  for  H^. 


The  condition  for  the  discontinuity  of  H at  y=0  is  obtained 

z 


by  combining  Eqs.  (101)  and  (117).  We  have 
2 

1 B fs^  cos  k 

X z 


1 2 

y (-1)  -r  (k  + — 5")  B (s)  cos  k x sin  k z cos  k b 


Z\  av  « 0%  J U 

2 2 mn 

m.n  c 

o 


V')  - ^ 2 

ktt 


r r (k' . k')  - (s^/cVk'^)e  (k'.k') 

I dk^  I dk^  { ^ ^ ^ — -—-M 


sp  k' 

o y 


-j(k^x  + k^z) 


(124) 
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where  E (k',  k')  and  E (k',  k')  can  be  determined  by  substituting 

Z X Z X X z 

Eqs.  (118)  and  (121)  Into  Eqs.  (115)  and  (104)  respectively. 

The  application  of  the  method  of  moments  to  Eq.  (124)  yields 
the  following  matrix  equation: 


[Z'(s)] 


a (s)  1 
pq  ' 


{•^(s)} 


(125) 


where  {<^(s)}  represents  the  excitation  column  matrix  which  Is  the  Inner 
product  of  ^j^(s)  and  the  testing  functions.  The  procedure  of  solution 
is  to  determine  {a  (s)}  and  (B  (s)}  from  Eq.  (125).  With  {a  (s)}  and 

pq  pq  pq 

(8  (s)},  (A  (s)}  and  {B  (s)}  can  be  found  from  Eqs.  (120)  and  (123) 

pq  mn  mn 

respectively,  which  are  then  used  to  obtain  the  field  distribution  Inside 
the  cavity  from  Eqs.  (98)-(101). 


VI-4.  Numerical  Techniques 

The  numerical  technique  outlined  In  Section  V for  determining 
{a(s)}  can  be  used  here  for  finding  {a  (s)}  and  (B  (s)}  from  Eq.  (125). 

pq  pq 

In  Section  V the  cavity  fields  were  expressed  In  terms  of  expansion  func- 
tions In  three  dimensions,  resulting  In  a [Z(s)]  with  large  dimensions, 
whereas  the  alternative  approach  developed  In  this  Section  expresses 
the  cavity  fields  In  relatively  few  cavity  modes  and  the  field  at  the 
aperture  Is  a two-dimensional  expansion.  Hence  the  matrix  [Z'(s)]  In 
Eq.  (125)  has  a much  smaller  dimension  than  [Z(s)]. 

In  applying  the  method  of  moments,  twenty  subsections  were 
used  for  the  aperture.  Rectangular  pulses  were  chosen  as  the  expan- 
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slon  function  in  the  x-dlrectlon  and  triangular  functions  as  that  In 
the  z-dlrectlon.  The  starting  points  of  the  root-searching  procedure 
for  det[Z'(s)]  ■ 0 were  the  resonant  frequencies  of  the  cavity  modes 
without  the  slot  aperture.  Unfortunately  [Z'(s)],  though  of  a smaller 
dimension,  was  also  Ill-conditioned  and  convergence  to  pole  singularities 
was  not  apparent  In  many  cases.  This  difficulty  could  not  be  alleviated 
by  changes  in  the  number  of  subsections  taken  at  the  aperture.  (This 
was  tried  without  success.)  An  effort  was  made  to  reduce  computation 
errors  by  increasing  the  magnitude  of  each  element  of  [Z'(s)]  by  a 
hundredfold.  This  yielded  the  two  layers  of  pole  singularities  shown 
In  Fig.  8.  The  convergence  was  very  slow  and  a time  limit  of  30  minutes 
was  set  for  locating  each  pole  on  the  Honeywell  6000  computer.  As  can 
be  seen  from  Fig.  8,  some  of  the  points  do  not  fall  on  the  smoothed 
layers,  indicating  computing  inaccuracies.  Much  time  was  spent  on  iso- 
lating possible  causes  for  the  inaccuracies  without  noticeable  success. 

The  computer  program  was  checked  and  rechecked  many  times.  Possible 
contributing  factors  to  the  slow  convergence  and  numerical  Inaccuracies 
are  round-off  errors,  internal  resonances  of  [Z'(s)],  and  existence  of 
singularities  other  than  simple  poles. 

The  normalized  natural  mode  for  the  first  pole  singularity 
of  the  first  layer,  Sj^^^  = (-0.121  + j2.039)c,  is  plotted  in  Fig.  9. 

The  natural  modes  for  the  other  pole  singularities  have  also  been  plotted. 
They  are  not  reproduced  here  because, except  for  slight  differences  in  the 
phase,  they  are  surprisingly  similar.  Explanations  for  this  phenomenon 
and  for  the  different  curvatures  of  the  loci  in  Fig.  8 as  compared  with 
those  in  Fig.  3 are  lacking. 
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Attempts  were  made  to  obtain  the  transient  electric  field  at 
typical  points  Inside  the  cavity  for  an  Incident  step  excitation,  but 
the  resulting  response  violated  causality  and  showed  peculiar  fluctu- 
ations and  therefore  does  not  appear  to  warrant  Inclusion  In  this 
report . 

As  a last  resort,  an  FFT  (fast-Fourler-transform)  algorithm 
was  applied  to  obtain  a time-domain  response  directly  from  the  frequency- 
domain  response. 

E(t)  = I E(s)e®*'ds  (126) 

C 

Substituting  s = jo),  we  have 

.00 

E(t)  = i I E(ja))e^“*^daj 

E(ja))e^“*^da)  (127) 

"■  -n 

Let  tj^  = kAt  = kT/N,  f^  = nAf  = n/T.  Equation  (127)  can  be  approximated 
as 

E(k)  = (Af)  I E(n)e^^’^"’^^^  (128) 

n=0 

where  E(k)  and  E(n)  stand  for  E(kAt)  and  £(j2TrnAf)  respectively.  In 

order  to  find  E(n),  Eq.  (12"5)  must  first  be  solved  for 

{0  (jw)}  which  are  used  for  obtaining  {A  (iio  )}  and  {B  (iw  )}  from 

pq  nr  mn  n mn  n 

Eqs.  (120)  and  (123).  The  frequency-domain  field  components  can  then 
be  determined  from  Eqs.  (98)-(101). 
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Two  FFT  runs  were  made  on  the  problem  with  perpendicular  polari- 
zation: one  for  evaluating  £(n)  at  257  frequency  points  and  one  at  513 
frequency  points.  Unfortunately  the  plots  for  E(k)  also  violated 
causality  and  showed  peculiar  fluctuations. 


VII.  CONCLUSION 


In  this  report  the  transient  excitation  of  a cavity  through  an 
aperture  is  used  to  illustrate  the  state-space  solution  of  transient 
electromagnetic  problems.  Depending  on  whether  the  E-  or  the  H-formulation 
is  used,  the  effect  of  the  aperture  is  accounted  for  by  a magnetic  current 
in  invoking  the  induction  theorem  or  by  an  electric  current  in  invoking  the 
equivalence  theorem.  The  governing ' second-order  differential  equation  is 
converted  into  a set  of  first-order  normalized  state  equations  by  defining 
three  new  state  variables  in  addition  to  the  appropriate  vector  potential. 
These  state  equations  subject  to  the  associated  boundary  conditions  are 
solved  by  the  method  of  moments.  The  cavity  region  is  first  divided  into 
subsections  and  the  fields  within  the  cavity  expressed  in  terms  of  appro- 
priate expansion  functions  with  time-dependent  coefficients.  The  fields 
in  the  half-space  outside  the  cavity  are  represented  as  superpositions  of 
plane  waves.  At  the  aperture,  the  cavity  and  external  fields  are  properly 
matched.  Inner  products  are  taken  with  testing  functions  and  the  first- 
order  equations  are  converted  into  matrix  equations  containing  the  expan- 
sion coefficients  as  unknown  column  vectors.  Evaluation  of  these  expansion- 
coefficient  vectors  leads  to  the  determination  of  field  distributions.  The 
procedure  for  evaluating  a typical  coefficient  vector  by  the  singularity- 
expansion  method  has  been  outlined. 

This  solution  procedure  was  applied  to  the  problem  of  electromagnetic 
excitation  of  a rectangular  cavity  through  a slot  aperture.  Two  cases  were 
considered:  an  E-formulatlon  for  parallel  polarization  and  an  H-formulatlon 
for  perpendicular  polarization.  For  the  case  of  parallel  polarization,  some 
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significant  pole  singularities  (natural  frequencies)  In  the  upper  left 
complex  plane  were  found  for  a sample  cavity.  The  corresponding  natural 
modes  were  determined  and  combined  to  give  the  electric  Intensities  as 
functions  of  time  at  two  locations  in  the  cavity  for  a step-excltatlon. 

Both  the  shape  and  the  magnitude  of  the  penetrated  electric  field  appear 
reasonable.  However,  It  would  not  be  wise  to  claim  absolute  accuracy 
because  of  the  Inherent  noise  and  instability  in  the  numerical  procedure 
and  because  of  a lack  of  a theoretical  proof  that  there  would  be  no  multi- 
ple poles  or  branch  points. 

For  the  case  of  perpendicular  polarization,  the  system  Impedance 
matrix  resulted  from  the  application  of  the  three-dimensional  moment 
method  vas  ill-conditioned  and  it  was  not  possible  to  determine  the  system 
singularities.  In  order  to  avoid  the  numerical  difficulties  inherent  with 
large  matrices,  an  alternative  analytical  approach  for  expanding  the  fields 
inside  the  cavity  in  terms  of  TM  and  TE  cavity  modes  with  unknown  coef- 
ficients was  developed.  This  approach  effectively  changes  the  three- 
dimensional  moment  method  into  a two-dimensional  one,  thus  reducing  the 
dimensions  of  the  resulting  matrix.  Some  pole  singularities  were  found. 
However,  the  system  matrix  was  also  unstable  and  convergence  to  the  singu- 
larities very  slow.  The  matrix  appears  to  contain  internal  resonances. 

An  FFT  algorithm  was  written  to  determine  the  transient  time-domain  response 
directly  from  the  frequency-domain  response.  Unfortunately  the  numerical 
results  violated  causality  and  showed  peculiar  fluctuations. 

As  mentioned  previously,  numerical  anomalies  appear  in  even  the 
simplest  cases  and  preclude  the  determination  of  the  natural  frequencies 
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of  well-defined  geometric  structures  [23], [24].  Numerical  difficulties 
in  cavity-backed  slot  aperture  problems  are  really  not  very  surprising. 
Previous  studies  [25],  [26]  on  the  penetration  of  transient  electromag- 
netic excitation  through  apertures  in  an  infinite  ground  screen  (without 
a cavity)  have  not  yielded  numerical  results.  Various  difficulties  were 
encountered  in  attempts  to  obtain  self-consistent  pole  locations  because 
of  numerical  instability  [26] . The  problem  is  vastly  more  complicated 
when  there  is  a cavity  behind  the  aperture.  It  has  been  shown  [27]  that 
the  structure  of  a conducting  cylinder  within  a parallel-plate  region 
gives  rise  to  two  types  of  singularities;  namely,  poles  and  branch  points, 
in  the  complex-frequency  plane.  Whether  a cavity-backed  slot  has  branch 
points  or  multiple  poles  is  an  unanswered  question.  It  is  felt  that  a 
concerted  effort  to  find  a definitive  answer  to  this  question  would  be 
a worthwhile  undertaking. 
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Fig.  6.  Response  to  Step  Excitation  (Parallel  Polarisation) 
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